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ABSTRACT 

We perform a detailed study of the electric and chromoelectric dipole co- 
efficients in i3 ^ Xg^ decay in a supersymmetric scheme with explicit CP 
violation. In our analysis, we adopt the minimal flavor violation scheme by 
taking into account the tan /3-enhanced large contributions beyond the lead- 
ing order. We show that the coefficients can deviate from the SM prediction 
significantly in both real and imaginary directions. Experimental bounds 
still allow for large deviations from the SM predictions for both dipole coef- 
ficients such that the CP asymmetry is as large as ±8%. There are further 
implications of these coefficients for the charmless hadronic and semileptonic 
B decays. As a direct application of our analysis, we have discussed A7 
decay. 







The Operator Product Expansion (OPE) combined with the heavy quark effective theory 
(HQET) forms the basic tool in analyzing the decays as well as productions of hadrons (See 
the review Basically, the effective Hamiltonian describing the scattering processes can be 
expanded in a series of local operators (whose hadronic matrix elements constitute the long- 
distance effects) with Wilsonian coefficients (which are generated by the short-distance effects) . 
Among all the hadronic scattering processes the rare ones are particularly important, as the 
contributions of the standard electroweak theory (SM) and those of the 'new physics (NP)' 
arise at the same loop level thus suffering no relative loop suppressions. Furthermore, decays of 
the 6-flavored hadrons {B, B*, A},, ■ ■ •), compared to strange and charmed ones, are especially 
important as for such systems the HQET is fully applicable, and via the OPE, one can both 
test the SM and search for possible NP effects, by confronting the associated Wilson coefficients 
with the experiment. 

In this work we will analyze the electric and chromoelectric dipole coefficients (denoted 
hereafter by C7 and Cg, respectively) describing the short-distance physics effects in rare B 
decays {e.g. B — > Xsj, B — > K*{892)j and B — > Xgt^t^ , • • •) 0], by taking into account existing 
experimental results. The electric dipole coefficient C7, rescaled to fi = rrib level, is directly 
constrained by the experimental result on B ^ Xs"f ||3|. However, the situation for the coefficient 
Cs is obscured by the fact that the gluonic decay 6 — > is not directly accesible in experiments. 
However, this very coefficient plays an important role in the charmless hadronic B decays. For 
example, the theoretical predictions for the inclusive semileptonic decay rate r(6 — > ceu) and the 
charm multiplicity in B meson decays are significantly higher than the experimental results (See 
1^] and references therein). The most plausable solution to this discrepancy stems from possible 
enhancement of the chromoelectric coefficient Cs by NP effects. Moreover, it is the relative phase 
between and C7 which determines the CP asymmetry in B ^ Xg'j, which will be measured 
in near-future B factories [^. Consequently, it is essential to determine the size and phase of 
the gluonic coefficient in regions of the parameter space of NP where the theoretical predictions 
for B — > X<j7 agrees with the experiment. 

In what follows, we will take low-energy minimal supersymmetry (SUSY) with explicit CP 
violation as the NP candidate. We will adopt the minimal flavor violation scenario (MFV), 
take into account the tan /3-enhanced large contributions beyond the leading order (LO). The 
inclusive mode B — > Xsj has already been analyzed within such a scheme by Q (in CP- 
conserving SUSY), and has been furthered by 0] to the CP-violating SUSY concluding a sizable 
CP asymmetry, which can compete with the experiment in near future. However, in both and 
10] the size and phase of Cg, its correlation with C7, and its interdependence with the branching 
ratio, as well as its effects on the CP asymmetry have not been reported in detail. In particular. 
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given the post-LEP bound on tan (3 3.5, it is necessary to have a detailed knowledge of Cg 
in this portion of the SUSY parameter space. The main goal of this work is to determine the 
size and phase of the chromoelectric coefficient Cg within the CP-violating SUSY at beyond- 
the~leading-order (BLO) precision in regions of the parameter space allowed by the existing 
B ^ Xsj constraint. 

The inclusive decay B Xsj is well approximated (within at most 10%) by the partonic 
decay b ^ which is described by the effective Hamiltonian 

Heff = -^V,*,Vt,j2C.{fi)0,{^) , (1) 

where V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and the operator basis 04=1,.. .,8 
is defined in fs], ^. This OPE for the Hamiltonian separates the long-distance (the matrix 
elements of the local operators Oi) and short-distance (associated Wilson coefficients Cj) at any 
scale /X G {rrib, Mw). Moreover, HQET approximates the inclusive rate by the partonic one (all 
terms being of the order 0{A.QCD/^b) and higher ones are negligible) @. 

Evolution of the Wilson coefficients from ^ = M\y down to /i = mi, level is governed by the 
standard QCD RGEs: 

C2(m,) = l/2(,?-12/23+^6/23) 

C7(m5) = C7(M^y)r?16/23+C8(MH.)^(??l^/2=^-??l^/23)^C^(^^)^;^.^r, 

1=1 

8 

C8(mfe) = Cs{MwW^'^^ +C2{Mw)Y.9iV'' (2) 

i=l 

where 77 = as{MyY) / as{mi,) , and the numerical coefficients /ij , rj and Qi are given in |p. 

The initial values for the QCD RGEs, C2Jfi{M^r), depend on details of the short-distance 
theory at ~ Myy. In standard electroweak theory for instance, one finds C2(m-fe) = 1.023, 
CriiT^b) = —0.312 and Cs{mh) = —0.148 at BLO precision Consequently, the NP effects 
can appear in various ways: (i) There can be observable deviations from these numbers with or 
without sign change, or (ii) the coefficients can take complex values. Each type of departure 
from the SM prediction implies certain aspects of the weak-scale NP effects. For instance, for 
(ii), it is obvious that the NP brings new sources of CP violation, and necessarily, the CP 
asymmetry of the decay deviates from the SM prediction (< 1%) |lO| . 

In what follows, we will take SUSY with explicit CP violation as the NP candidate, and 
concentrate on the results of 0] where the Wilson coefficients were computed at NLO precision 
for those threshold effects enhanced at large tan/3. Within this framework C2{My[r) = 1 as in 
the SM, but the two dipole coefficients Cj^si^w) are significantly modified compared to the LO 
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results With the MFV scheme, only chargino-top squark, charged Higgs-top quark and 

H^-boson-top quark loops give significant contributions 

C7,8(Mm/) = CY^siMw) + C^,s{Mw) + ClsiMw) (3) 

In CP-violating SUSY, at LO precison C^^(Miy) are always real; they do not contribute to 
CP-violating observables, such as the CP-asymmetry, in the decay. However, the chargino 
contribution is complex due to the fi parameter (having finite phase and the stop trilinear 
coupling At (having finite phase (Pa) with approximate structure e*^'^''^'^'*^ [11|. 



With BLO precision, however, there are finite threshold corrections to each piece in (^) such 
that all three contributions C^q{Mw), Cj^g^{Mw) and C^giMw) are now complex. Moreover, 
larger the tan /? larger their imaginary parts, so that even naively one expects CP-violating 
effects to be enhanced at large tan/3. Indeed, as reported in Q, the BLO CP asymmetry 
is significantly larger than the LO one at sufficiently large tan/3. Although the asymmetry 
remains ^ 8% in both cases there occurs an enhanced sensitivity to tan /? for the BLO case. 

In the following we will perform a numerical study of the electric and chromoelectric dipole 
coefficients, and discuss their phenomenological implications. In the numerical analysis, we take: 
(i) the light stop t2 and the charged Higgs are degenerate and weigh close to the weak scale, 
= = 250 GeV; (ii) the sfermions of first two generations are heavy enough, so that 
one can neglect their contribution to ^ Xsj (this is a viable way of suppressing the one-loop 
EDMs §); (m) the SU(2) gaugino, the right-handed sbottom and the heavy stop are heavy, 
M2 = M~ = fht,^ = 1 TeV, and form the SUSY breaking scale; (iv) the stop and sbottom 
trilinear couplings have the same phase, 9a,, = 9 At, and the latter is degenerate with the fi 
parameter = \Af)\ = 150 GeV (sbottom parameters are needed for 2-loop EDM calculations 
||l3|| ); and finally (v) the light stop is dominantly right-handed to agree with the electroweak 
precision data, hence the stop mixing is to be sufficiently small (9^ = 7r/20). Moreover, we 
vary tan (3 from 10 to 50 and the phase 4'A,fi from to vr in forming the scatter plots. The 
parameter space mentioned here has been determined after trying several combinations, and it 
is one those points yielding large CP-violation in the system. One particularly notices that the 
lightest chargino is a Higgsino so that CP-violation via chargino-stop contribution is enhanced. 
In minimal supergravity, for instance, the lightest chargino is SU(2) gaugino and thus CP- 
violation is very much suppressed. 

Depicted in Fig. 1 is the tan /? dependence of Re[C7(m;,)] and Im[C7(m;,)] for 10 < tan/3 < 50 
and < (j)fj_,(j)A ^ TT. When tan/3 ^ 10, Re[C7(m;,)] is close to the SM value and Im[C7(m5)] 
is evenly distributed around the origin, being consistent with zero. However, as tan [3 rises 
towards larger values, so does Re[C7(mf,)], crossing zero around tan/3 ~ 27. When Ian (5 ^ 27, 



3 




10 15 20 25 30 35 40 45 50 
tan P 




10 15 20 25 30 35 40 45 50 
tan /3 

Figure 1: The tan/3 dependence of Re[C7(m;,)] (upper window) and Im[C7(m6)] (lower window), 
for the values of 0^ and 4)a varying from to tt. 

^e[Cj{rai,)] takes both negative and positive values, and it can be as large as 0.4. The imaginary 
part of Cy, however, rises with increasing tan/3 in absolute magnitude for both negative and 
positive directions in an approximately symmetric manner. For tan (3 ^ 35, its distribution is 
levelled with a value swinging between -0.35 and +0.35 smoothly. Clearly, for large enough tan/3 
there are regions of the parameter space where Im[C7(m;,)] is not consistent with zero, signalling 
therefore a clear sign of the NP effects. An approximate expression valid for large tan (3 can 
be given as Cj ~ 0.4 e^'""/*, which is far away from the SM prediction in both size and phase. 
However, there are regions of the parameter space where is pure imaginary, is pure real or 
vanishes exactly. 

Similar to observations made for Fig. 1, one can discuss the Wilson coefficient Cg using Fig. 2 
where its real (upper window) and imaginary (lower window) parts are separately plotted against 
tan (3 when (/>^,a vary from to vr. One notices that, unlike Cj, for low tan (3, Re[C8(?Tife)] deviates 
from its SM value. This stems from the fact that Csinib) is directly proportional to CsiMw) (up 
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Figure 2: The tan j3 dependence of Re[C8(m6)] (upper window) and Im[C8(mft)] (lower window), 
for the values of (j)^ and (j)A varying from to vr. 

to small corrections proportional to C2{My\r))^ and therefore the NP effects at short distances 
are directly reflected to the hadronic scale. As tan/? rises to larger values, both Re[C8(mfe)] 
and Im[C8(mb)] gradually increase where the former crosses zero around tan/? ~ 35. In large 
tan (3 regime, Re[C8(mfe)] can take positive values only in a small portion of the parameter space 
whereas \m\C^{rrn,)] swings between -0.1 and 0.1 almost evenly. One notices that, as in C7, in 
certain regions of the parameter space C^{Mw) can be pure real, pure imaginary or just vanish. 
Those points where both C^ and Cs vanish are particularly interesting, as in this case one has 
to saturate the experimental bounds on B ^ X^^ via the chirality-flipped Wilson coefficients, 
implying large contributions due to the gluino exchange |jl^. As we are working in the MFV 
scheme, such effects are obviously beyond the scope of our discussion. 

Depicted in Fig. 3 is the dependence of Im[C8(m{,)] on Re[C7(m;,)]. Here the size and shape of 
the ellipsoidal region is determined by the accuracy of the experimental results on BR(i? — > ^^7) . 
The main NP effect occurs in pushing ^e[C-j{m},)] to larger values compared to the SM prediction 
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instead of smaller ones. The small region around Re[C7(mb)] ~ —0.3 corresponds to the SM 
validity domain where imaginary parts of both coefficients remain around zero. 

0.15 
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Figure 3: Dependence of Im[C8("ife)] on Re[C7(mfc)]. The width of the ellipsoid corresponds to 
present experimental accuracy for BR(i? ^sT) • 

Fig. 4 illustrates the dependence of A(jp{B Xgj) on Re[Cs{mh)] (upper window), and 
on Im[C8(m;,)] (lower window). As we can see from the upper window, the asymmetry takes 
the largest value ±8% when Re[C8(?7ife)] is in the range of the SM prediction for Csimb). One 
notes that in a more general model CP asymmetry can be as large as 20% [l^]. Moreover, in 
certain small corners of the SUSY parameter space the CP asymmetry can be enhanced by 
a factor of 2 |jl^. When Re[C8(mb)] takes larger or smaller values than the SM prediction, 
the asymmetry gradually drops to the corresponding SM prediction ~ 1%. The dependence of 
the CP asymmetry on Im[C8(m6)], however, shows that the maximal values are attained when 
Im[C8(mfe)] ~ ±0.1, which is far away from the SM predicton. In conclusion, the asymmetry is 
roughly an order of magnitude larger than what is expected in the SM, and this happens when 
Re[C8("T'f))] remains close to the SM prediction, and Im[C8(r7ib)] is large and has a similar size as 
the real part. 

A closer comparative look at the figures suggests that the CP asymmetry, which will be 
measured in near-future B factories with increasing precision, is maximal when |C8(7Tifc)| ~ 



|C7(mfe)|. This results confirms earlier predictions |16| where it was already shown that Aqp{B 
X,7)~10% \Cs{mb)/C7imb)\. 

Even after including the NLL corrections to the semileptonic inclusive B decay B Xev, 
the theoretical predictons for the branching ratio and the charm multiplicity turn out to be larger 
than the experimental result |17, ^. Therefore, it is conceivable that possible enhancements in 
Cs can account for the existing discrepancy between the theory and the experiment. As the 
numerical analyses above show, it is possible to significantly shift this coefficient in both real 
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Figure 4: The dependence of Aqp{B X^j) on Re[Cs{mi,)] (upper window) and on Im[C8(mfe)] 
(lower window). 

and imaginary directions compared to the SM prediction. Therefore, the experimental result 
on i? — > Xg'j allows for large deviations in chromoelectric as well as electric coefficients at 
large values of tan /? where the SUSY threshold corrections are important. Concerning the 
semileptonic B decays, the exclusive channel B — > Xsi^£~ is another example having both 
theoretical and experimental importance. This follows from the fact that the forward-backward 
asymmetry Apg of this decay vanishes at a specific value of the dilepton invariant mass [|l8| in 
a hadronically clean way. This zero of Apg occurs at the point mf^ ~ — m^Re [Cj{mi,)/CQ{mb)] 
where Cg{mij) is the four-fermion operator coefficient (not computed here). The important 
point is that the position of the zero shifts in accord with the NP contributions to the Wilson 
coefficients. Suppose that C-/{mh) alone is shifted in complex direction, then it is clear that the 
critical value of the dilepton mass mj^ is shifted back and forth, depending on the parameter 
values @. 

There are other processes where the radiatively corrected Wilson coefficients play an im- 
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portant role. As an application, for instance, we analyze the bottom baryon radiative decay, 
Afe — > A7, which is again dominated by the mechanism b — > 57. We will particularly concentrate 
on the dependence of the braching ratio on the Wilson coefficients at the weak scale. 
The decay rate of A;, — > A7 has been computed in , and its expression is given by: 




r(A,^A7) = |a|^ + |6|^ , (4) 



with 



Gf e 
72 8^ 



2Cr{mb)mbVM { /f ^^(0) - /^''(O) 



(5) 



The form factors in a and b coefficients, depend on the momentum transfer as follows jl^, pO| : 




/my j V ~ H l™'A 



where = (jnk^ — vn^) and, q = — pA, with n=l and n=2 representing the monopole 
and dipole contributions. Here, my and ruA are the pole masses of the vector and axial vector 
mesons, respectively. 

In the numerical analysis, following |19, 2C], we let fi^''{Qm) — 9i^*'{lm) — 0.64, 
92^''{lrn) ~ —0.10, and f2^^{Q^) = —0.31 for the values of the heavy light form factors. Addi- 
tionaly, we take 'mi,{mi,) = 4.25 GeV, niA^ = 5.624 GeV, my = 5.42 GeV, m^ = 5.86 GeV, and 
r(Ab) = 1.23 X 10~^^sec. In forming the scatter plots, we vary tan (3 from 10 to 50 and the phase 
4'A,fj. from to vr. 

Depicted in Fig. 5 is the dependence of the branching ratio of A^ A7 (^(Af, A7)) on 
Re[Cs{M\Y)] (upper window) and Im[C7(Miy)] (lower window), for 10 < tan/3 < 50 and < 

4'A ^ where only the monopole dependence for baryon form factors is considered (n=l). 
As is seen from both windows, Re[C8(Mvy)] takes negative values up to ^(A^ — > A7) ~ 1.7x 10^^, 
whereas Im[C8(Mvy)] is evenly distributed around the origin in this interval. Re[C8(Mvy)] starts 
to take both positive and negative values, when I3{Af, A7) > 1.7 x 10~^. It can take only 
positive values for 13{Aiy A7) ^ 2.7 x 10~^, and can be as large as ~ 0.2. On the other hand, 
for 13{Ai, A7) > 1.7 X 10^^, Im[C8(Mvi/)] increases in both positive and negative directions in 
a symmetric manner. Clearly, there are certain regions of the parameter space where Cs{M]y) 
can be pure real, pure imaginary as well as it just vanishes. 

In the recent work |19], it has been shown that B{Ah — > A7) has a magnitude of 1.9 x 10^^, 
when only the monopole q'^ dependence of the baryon form factors is considered (n=l). This 
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1 1.5 2 2.5 3 3.5 

S(Aj ^ A7) X 10^ 




1 1.5 2 2.5 3 3.5 

B(At A7) X 10^ 

Figure 5: The dependence of B{Ai, A7) on Re[C8(Mvi/)] (upper window) and on Im[C8(Mvy)] 
(lower window), for the values of (j)^ and varying from to vr. 

decay has also been considered in ||2l| with the predicted branching ratios in the range of 
(1.2 — 1.9) X 10~^. As one can see from both windows of Fig. 5, ^(A;, — > A7) varies in the range of 
(1— 3.2)xl0~^, and the maximal values of i3(A;, — > A7) are attained when Re [C8(Mi4/)] ~ 0.2, and 
Im[C8(Mvy)] ~ ±0.2, which are away from the S.M prediction for Cs{M]y). A closer comparative 
look at the figure suggests that when Cs{Mw) is in the range of the S.M prediction, for instance 
CsiMw) = -0.086 ||], ^(Af, A7) drops nearly to ~ 1.9 x 10"^, which is consistent with the 
results of S.M[^]. However, one notes that at this particular value of the branching ratio, there 
are other solutions in the parameter space where CsiMw) can deviate from its S.M value not 
only in sign, but also in magnitude, as it can be complex. 

We would like to note that the analysis of |jl9| has been carried out in the context of the S.M, 
and C7(mfc) has a fixed value {Ci{mb) = —0.312). However, in our work, C7 and Cs are complex, 
and as we have shown in Figs. 1 and 2, they vary in a wider range. For instance, approximate 
expressions of 6*7(8) valid for large tan/3 can be given as C7 ~ 0.4 e^^'^^^, and Cg ~ 0.1 e^*'^/^. 
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which are far away from the S.M predictions in both size and phase. Therefore, having a wide 
range of parameter space, the model under concern gives new allowed regions. It is natural that 
there are regions in the parameter space where our prediction for ;B(Af, — > A7) agrees with that 
of S.M. 




1 1.5 2 2.5 3 3.5 

B(Aj ^ A7) X lO''^ 




1 1.5 2 2.5 3 3.5 

B(At A7) X 10-5 

Figure 6: The dependence of i3(Afe — > A7) on Re[C7(MvK)] (upper window) and on Im[C7(Mi4/)] 
(lower window) , for the values of (pfj, and varying from to vr. 

Similar to the observations made for Fig. 5, one can discuss the Wilson coefficient C7 using 
Fig. 6, where its real (upper window) and imaginary (lower window) parts are separately plotted 
against B{Ah A7), when (/>^,a vary from to tt, and only monopole contribution is considered 
(n=l). As B{Ah A7) rises to larger values, both Re[C7(Mvy)] and Im[C7(Mvi/)] gradually 
increase, where the former crosses zero around B{Ab A7) ~ 1.3 x 10"^. Re[C7{Mw)] starts to 
take positive values only when B{Ah — > A7) > 2.4 x 10~^, whereas Im[C7(Mty)] swings between 
±0.6. Like Cg, there are observable deviations from the S.M prediction with or without sign 
change, and it can take complex values as well. However, a comparative glance at the figure 
shows that, when Ct{Mw) is in the range of the SM prediction, for instance C7{M\y) = — 0.161[^, 
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1 1.5 2 2.5 3 3.5 4 

B(Aj ^ A7) X 10*5 

Figure 7: The dependence of B{Ah — > A7) on Aqp{B Xg'y) for n=l (upper window) and for 
n=2 (lower window) 

Depicted in Fig. 7 is the dependence of B{Af, A7) on A(jpiB ^sj), when only the 
monopole (upper window) and the dipole (lower window) dependences are considered. As is 
noticed from the upper and lower windows S(Af, — > A7) is of the orders of 10~^ and 10^^, for 
the monopole (n=l) and dipole (n=2) contributions, respectively. In both cases ^(Af, A7) 
behave similarly. However, it decreases by an order of magnitude in the dipole case, and the 
maximal values of Aqp{B Xg'y) are obtained when B{Ai, A7) ~ 1.5 x 10~^(1.9 x 10~^) for 
n=l (2). In the recent work of ||l^, it has been shown that B{Ai, A7) = 2.3 x 10~^ for n=2. 
As we can see from the upper and lower windows, the asymmetry nearly takes the largest value, 
when B{Ab A7) is in the range of the SM prediction |jl|, |2|, ||]. When ^(A;, A7) takes 
larger values than the SM prediction, the asymmetry gradually drops to the corresponding SM 
prediction ~ 1%. 

In conclusion, we have computed the dipole coefficients Cj^s in SUSY with explicit CP viola- 
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tion with special emphasis on large values of tan (j. We have shown that the present experimental 
bounds oYiB ^ Xg"f allows for large deviations in the Wilson coefficients (with respect to the SM 
prediction) in both real and imaginary directions. The CP asymmetry in the decay is enhanced 
by an order of magnitude, thanks to especially the SUSY threshold corrections. The allowed 
deviations from the SM values can account for (being a plausable hypothesis) the discrepancy 
between the experiment and theory for the semileptonic B decays. As an illustration, we have 
discussed Af, — > A7 decay. 

M. B would like to thank the Scientific and Technical Research Council of Turkey (TUBITAK) 
for partial support under the project, No:TBAG2002(100T108). 
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